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1 Introduction 



In recent years there has been tremendous interest in the spin models of one-dimensional magnets because materials 
sciences have made a remarkable progress and a number of quasi-one-dimensional magnets have become available. 
On the one hand, the one-dimensional magnets may exhibit a wide variety of peculiar properties (e.g., gapless or 
gapped low-energy excitations, magnetization plateaus, spin-Peierls phases etc.). On the other hand, for the one- 
dimensional spin models specific theoretical tools can be applied (Jordan-Wigner fermionization, bosonization, Bethe 
ansatz etc.) and therefore a more comprehensive analysis can be performed. Although usually the Heisenberg model 
with modification is used to describe the properties of low-dimensional magnetic compounds, some of their generic 
features can be illustrated within the simpler framework of the XX model in a transverse field. 

In the present paper we study the effects of regular modulation of the exchange interactions on the magnetization 
and static susceptibility examining for this purpose the local magnetic properties of the spin-i XX chain in a transverse 
field. A model with the alternating exchange interactions arises while describing superlattices or annealed "bond 
impurities" || . Recent study of the alternating quantum (Heisenberg) spin chain with experimental applications has 
been reported in Quantum spin chains with regularly modulated exchange interactions attract much attention 
because of the magnetization plateaus which can be realized in such chains. Oshikawa, Yamanaka and Affleck Q 
considered a general quantum spin chain exhibiting rotational symmetry with respect to the direction of the applied 
uniform field and proposed the necessary condition for magnetization plateaus as p{s — m) = integer, where p is the 
period of chain, s is the spin value, m is the zero temperature magnetization per site. Much work has been done [5-9] to 
investigate magnetization plateaus in various quantum spin chains using approximate analytical approaches (mainly, 
bosonization plus renormalization group analysis or strong-coupling limit plus perturbation expansions) and numerical 
ones (mainly, exact diagonalization or density-matrix renormalization group techniques). On the other hand, for a 
simpler spin-i XX chain the magnetization profiles can be calculated exactly for any finite period of nonuniformity 
p [|f0). Moreover, as we shall show below, it is also possible to examine rigorously a nonuniversal behavior of the 
local magnetizations of such chains. In addition, we can calculate the local static susceptibilities. The local magnetic 
properties of quantum spin chains can be probed experimentally. Thus, the precise analysis of the NMR line shapes 
can yield the on-site magnetizations in the corresponding quasi-one-dimensional compounds |pd| , |l2f . 

The paper is organized as follows. First we explain how the local magnetic properties of the regularly nonuniform 
spin-i XX chain in a transverse field can be calculated exactly. Then we analyze the obtained results for the zero 
temperature magnetization profiles from a viewpoint of the strong-coupling limit. We show to what extent the exact 
zero temperature magnetization profiles can be reproduced within the framework of the strong-coupling approximation. 
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Further we discuss the temperature dependence of local static susceptibilities. We end up with a brief summary. 



2 Continued fractions and the local magnetic properties 

Consider a nonuniform spin-i XX chain in a transverse field defined by the Hamiltonian 

N N 

H = J2 Vnsl + 2 £ I n + s« n sl +1 ) (1) 

7i—l n — 1 

N ( i\ N 

= ( s n s n ~ 7} ) + ^« ( s n s n+l + s n s n+l) ■ 



n— 1 v 7 n— 1 

By applying the Jordan- Wigner transformation one comes to spinless fermions with the Hamiltonian 

N / , x N 



H = 0„ ( c+c„ - - J + ^ I n (c+c n+ i - c„c+ +1 ) . (2) 



2, 

n— 1 v 7 n— 1 



We shall be concerned with the local on-site magnetization m n = (s n ) = (s+s n ) — | = (c+c„) — i and the local 
on-site static susceptibility Xn — 9 qq ■ Here the angular brackets denote the canonical thermodynamic average. 

Let us introduce the Green functions (l3) G% m {t) = T'^{±t){{c n (t), c+ }), Gf lm (t) = ± dwe-^G^u ± ie), 
e — > +0, and note that for the model defined by Eq. (||) the following continued fraction representation for GJ„ = 
G^ n (uj ± ie) holds 

w ± ie — S2„ — A„ - A„ 



A: 



A, 



J n-1 



w ± ie — r2„_i 



w±le-f!„. 

/ 2 



lo ± ie — fi n +i 



aj±le — £1^+2- 



For any finite period of varying of the Hamiltonian parameters Q n and I n the continued fractions in Eq. (^) become 
periodic and can be evaluated by solving quadratic equations. As a result one gets the exact expressions for G^ n , 
(c+c„) = do; -g^ff , and hence for m„ and x„. Introducing the "local" density of states p n (u>) — ^^ImG^ 

the formulas for m n and % n can be rewritten as follows 

m n = -- / dwp„(w)tanh— , (4) 

* J— do * 

Here /3 = ^ is the reciprocal temperature. Knowing the local quantities p n (u)), m„ and Xn one immediately gets the 
(total) density of states p(ui) — Yl n =i Pn{u) an d the (total) magnetization and static susceptibility m = -4 Xm=i m « 
and X = E n =i Xn, respectively. 



Following the procedure outlined above for the calculation of p n (u>) (some further details can be found in |T(J 
where, however, only the total density of states p(u>) was considered) one can easily derive p n {uj) for chains provided 
the period is not too long. For example, for a chain of period 3 one finds 

j iJM, if ch>o, 

PnH = I ^ (6) 

0, otherwise, 

y n {i0) = (UJ - O n+1 ) (W - ^«+2) - In+l, 



c{uj) = 4i 1 2 ifif - (i 2 (w - n 3 ) + i| (w - Oi) + if (w - o 2 ) - (w - no (w - n 2 ) (w - n 3 )) 
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= -JJ( W_C i)' 

where Cj are the six roots of the equation C(u>) — 0. For a chain of period 4 the corresponding result is 

ir^nM if v( u )>o, 

,„U-. ^^M' (7 ) 
0, otherwise, 

W n = (w - O n+1 ) (u - O n+2 ) (w - O n+3 ) - I 2 +1 (w - - In+2 ( w - ^n+l) , 

v{uo) = uli^ilil - {{uj - no (w - n 2 ) - n 3 ) (w - o 4 ) 
-i 2 (w - n 3 ) (w - n 4 ) - i 2 2 (w - Hi) [u - n 4 ) 
-if (w - no (w - n 2 ) - 1 4 2 (w - n 2 ) (w - n 3 ) 

+I 1 2 lf + l|l|) 2 = -n(a;-d,) ) 

i=i 

where dj are the eight roots of the equation T>(u>) = 0. The analytical calculation of p n (ui) for chains of larger periods 
can be easily implemented on a small computer (see the results for periods 6 and 12 presented below). 



3 Magnetization: exact results 

Let us discuss the local magnetic properties of the spin model ([!]) induced by regular nonuniformity. In what follows we 
restrict ourselves to a case of the uniform field n„ = n and assume that I„ = 7(1 + 8 n ) where S n is taken either in the 
form 5 n = 5 (6 rh i - 5 na + <5„,i +p - $ n ,2+ P + <W+2 P - 8 n>2+ 2 P + • • • ) or in the form 5 n = -5 cos ^p. The parameter S 
controls the deviation from uniformity and the parameter p is the period of modulation of the exchange interactions. 
The modulation of exchange interactions may be interpreted as a result of the lattice distortion. For example, the 
first ansatz for S„ corresponds to a displacement of the second site towards the first one, of the (p + 2)th site towards 
the (p+ l)th one and so on. The second ansatz corresponds to the displacements of sites giving by cosine with period 
p. Some relevant lattice configurations for p = 3, 4, 6, 12 are shown in Fig. 1. In our illustrations presented below we 
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Figure 1: Lattice distortions of the form <5„ = 8 (6 n ,i — 8 n ,2 + S n ,i+ P — S„ t 2+ P + S n ,i+2 P — 8„,2+2 P + ■ ■ ■ ) (a, c, e, g) and of the 
form 5 n = -8 cos ^ (b, d, f, h) for p = 3 (a, b), p = 4 (c, d), p = 6 (e, f) and p = 12 (g, h). 

put for concreteness I — 1, 6 — 0.6 (except Figs. 4, 5 where S varies from 0.3 to 0.999). 

In Fig. 2 the zero temperature dependences of local magnetizations m n on field Cl are depicted. Similarly to the 
total magnetization m, the local magnetizations m n in regularly alternating XX chains exhibit a step-like dependence 
on the applied field. From the mathematical point of view this is the consequence of a splitting of the initial fcrmion 
band (p n (uj) = i— ======= if fi — 2|7| < oj < ft + 2\I\ and p n (uj) — otherwise) into several subbands, the edges 

of which are determined by the roots of equations C(ui) — (||), T>(ui) = (0), etc.. As a result the zero temperature 
dependence m n vs. fl [as well as m vs. Q] must be composed of sharply increasing parts (as can be seen from Eq. (^) 
they appear when to = remains inside a subband while f2 increases) separated by the horizontal parts (they appear 
when lu = remains outside a subband while Q increases). The characteristic fields at which the horizontal parts of 
the dependences m„ vs. £1 begin and end up (and therefore the plateau lengths) are the same for all m„ and coincide 
with the values of such fields for m, since in all the cases those fields are determined by the roots of equations C(u>) = 0, 
T>{u>) = 0, etc.. However, the heights of plateaus m n are essentially site-dependent quantities. Moreover, the values of 
m n depend on 5. This is in contrast to the behavior of to, since the possible values of to are universal and do not depend 
on details of the intersite exchange interactions. As can be seen, for example, in Fig. 2b for the lattice shown in Fig. 
lb the magnetization at site 2 in weak and in moderate fields is always smaller than the magnetizations at sites 1 and 3 
and its value depend on S (compare the results for S = 0.6 seen in Figs. 2a, 2b, 2e and for 8 = 1 presented below). On 
the other hand, the sequence m ni < m n2 < . . . < to„ p may depend on the value of the applied field f2 as can be seen in 
Fig. 2e (and Figs. 2g, 2h) . For example, for the chain shown in Fig. le —mi < — TO4 < — TO2 < — tos < — me < — to 3 
for 0.7037 . . . < ft < 1.3578 ... but -mi < -m 2 < -m 6 < -m 5 < -m 4 < -m 3 for 1.5884 ... < ft < 2.0139 .... 

The mentioned property is also visible in Fig. 3 where the zero temperature magnetization profiles along chains for 
various fields are displayed. Moreover, the sequence m ni < m n2 < . . . < m Ilp for a certain Cl may depend on 5 as 
can be seen in Fig. 4. The results presented in Fig. 2 (3) indicate a relation between the periodically modulated 




Figure 2: Local magnetizations m n vs. field Q at zero temperature f3 = oo for the chains shown in Fig. 1 (/ = 1, 5 = 0.6). 
Panel a corresponds to the chain shown in Fig. la, panel b corresponds to the chain shown in Fig. lb, and so on; the curve 
denoted by n presents the local magnetization m n ; the bold curves present m. 
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Figure 3: Zero-temperature magnetization profiles {m n } at f2 = 0.5, 1, 1.5, 2, 2.5 (from bottom to top) for the chains shown in 
Fig. 1 (7 = 1, 5 = 0.6). Dashed lines indicate m„ — 0. Panel a corresponds to the chain shown in Fig. la, panel b corresponds 
to the chain shown in Fig. lb, and so on. 
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Figure 4: Local magnetizations m n vs. field ft at zero temperature for the chain shown in Fig. le for S = 0.3 (a), 0.6 (b), 0.9 
(c), 0.999 (d) (I = 1). 
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exchange interactions [or lattice distortion] and magnetization profiles at different sites (magnetization profiles along 
a chain for various fields), which for the considered spin-i XX chain in a transverse field can be traced rigorously. 
Such correspondence may be of experimental importance. Thus, the NMR spectra yield the distribution of the local 
magnetization (see, for example, JH], ^]). In particular, such analysis was performed to confirm that the lattice 
modulation in the incommensurate phase in CuGeC>3 has the form of a soliton lattice fllTf . 

4 Magnetization: strong— coupling approximation 

Although we are able to calculate the magnetization profiles for the periodic nonuniform spin-i XX chain exactly 
it is worth to consider those results from a viewpoint of the strong-coupling approximation, which was exploited in 
a number of papers devoted to spin chains j^] and ladders j^] . On the one hand, by turning to the strong-coupling 
limit, we can understand better the magnetization processes at zero temperature j3 — oo having a physical picture in 
terms of spins rather than in terms of Jordan-Wigner fermions. On the other hand, we can use the exact results as 
a testing ground to see how the approximate approach works. In the strong-coupling limit we put 5 = 1. In such a 
limiting case the chain of period p splits into noninteracting clusters each consisting of p sites. The zero temperature 
local magnetization in such limit is given by m n = (GS\s^\GS) where \GS) is the ground state eigenvector of the 
cluster Hamiltonian. The magnetization plateaus arise due to a change of the ground state with varying of the field. 

Consider, for example, the chain of p = 3 shown in Fig. fa. The relevant 3-site cluster Hamiltonian has the 
following eigenvectors and eigenvalues 

|1) = -!=| Ul 2 h) + litals) - Trials), |2) = -L=| T1T2I3) + T1I2T3) - Uhh), 

|3) = I lilals), |4) = -^=| IU2T3) - -)=\ hhla), 
|5) = T1T2I3) - T1I2T3), 1 6) - I T1T2T3), 

1 7) = --L| UUh) + I1T2I3) + Tiiais), |8) = -j=| T1T2I3) + TiiaTs) + ^1 liTaTs), (8) 

E 1 = ~n-VEi, E 2 = ^n-V5i, s 3 = -|o, E < = ~n, 

E 5 = ~n, E 6 = ~n, E 7 = -^n + V5i, E s = ^n + V5i, (9) 

respectively. As follows from Eq. (Q) for < Q < \fhl the ground state is |1) and therefore according to Eq. (^) 
mi = 0, 77i2 = —jk, fn?, = — § , and m = — J, that agreeing with the results for 5 = 0.6 plotted in Fig. 2a. If SI exceeds 
then the ground state is |3) and according to Eq. (|J) one gets mi = 777,2 = 7713 = m = — \. Similar arguments 
applied to the chain shown in Fig. Ib yield mi = 7773 = — 4, 7772 = if < ft < -^=1 and mi = 7722 = ma = — ^ if 



Let us consider further the chain of p = 6 shown in Fig. le. The eigenvalues and the eigenvectors of the relevant 
6-site cluster Hamiltonian can be easily found numerically. Thus, the eigenvector \S Z = —1), i.e. with the expectation 
value of S z = sf + . . . + Sg equals to — 1, which yields the lowest energy in this sector of S z , is 

\S Z = -1) = -0.16| liUaUUU) + 0.21| lihlaUUU) - 0.14| lihhUUU) 
+0.23| iihhUUU) -0.23| TiWsUlsTe) -0.07| UUzUhU) 
+0.061 lilatsUtsle) - 0.38| litalsUTsle) + 0.41| TilalsUTsle) 
-0.021 |i|2t3T4U 6 )+0.39| I1T2I3T45W -0.44| hhlaUhU) 

-o.25| UhhUUU) + o.28| hhhUUU) - o.06| T1T24.3i4J.5i6>, (10) 

whereas the eigenvector \S Z = —2), i.e. with the expectation value of S z equals to —2, which yields the lowest energy 
for given 5 Z , is 

\S Z = -2) = 0.38| liU.UUU) - 0.22| UlilsUUU) + 0.12\ UhlsUUU) 

-0.05| iiiatsiiieie) + 0.58| TiTsisUisie) - 0.67| TiiaisiUsie) (H) 

(only two digits after decimal point are preserved in the expansion coefficients in Eqs. ([!(]), ([H])). Calculating the 
expectation values of s z n with the help of Eq. (|Io| ) and Eq. (|ll| ) one immediately finds that (S z = — l\s^\S z = —1) > 
(S z = -l\s z 6 \S z = -1) and (S z = -2|s||5 2 = -2) < (S z = -2|s||5 2 = -2) just what was observed in Fig. 2e, 
however, for finite S = 0.6. A validity of the strong-coupling limit predictions (which are exact at S = 1) can be 
estimated from Fig. 4 where we plot the exact zero temperature magnetization profiles for the chain shown in Fig. le 
with S varying from 0.3 to 0.999. 

The results of the strong-coupling limit may be used to obtain the approximate zero temperature magnetization 
profiles when 5 is slightly less than 1. Consider, for example, the chain of p = 3 shown in Fig la. The Hamiltonian 
can be naturally split into a sum of 3-site cluster Hamiltonians (main part) and a part describing the inter-cluster 
interaction (perturbation) . Only the two lowest levels of the 3-site cluster Hamiltonian should be taken into account 
to describe the zero temperature magnetization profiles in the region between m = — \ and m = — h. We assume 
O > and hence the lowest relevant levels are |1) and |3). Introducing spin | operators a a attached to each cluster 
which act as a z |1) = — g|l), cr^ |3) = ||3) etc. one can find that the zero temperature magnetization profiles can be 
calculated with the help of the Hamiltonian of the uniform spin-i XX chain in a transverse field 

1=1 \ L V / J l + (l + d) l=l 
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Moreover, the spin operators at sites s z n are expressed through a z as follows 

s\ = — - — -er 2 
1 4 2 

2+(l + <5) 2 (l + Sf _ ff2 



4(l + (l + 5) 2 ) 2 (l + (1 + 5) 2 ) 



l + 2(l + £) 2 1 n „s 

4(l + (l + 5) 2 J 2(l + (l + ,5) 2 J 

Dotted curves in Fig. 5 show how this effective Hamiltonian works while 8 deviates from 1. We can compare the 
exact and approximate results. As S decreases the strong-coupling approximation predictions for the values of m n , 
characteristic fields, and detailed shape of profiles start to differ noticeably from the exact results. However, even for 
(5 = 0.6 (Fig. 5b) the strong-coupling approximation yields a reasonably good quantitative picture of the magnetization 
process in the region between m — — i and m = — ^. Evidently, for small 5 (Fig. 5a) the strong-coupling approach 
becomes worse and fails in the uniform limit S — ► as can be seen from Eq. (|l3|). Obviously, the strong-coupling 
approach can be used for other chains with regularly modulated exchange interactions (e.g. the Heisenberg chain), 
for which there are no exact results and the results presented in Fig. 5 may be of use to estimate the accuracy of the 
strong-coupling approximation. 



5 Susceptibility 

Last let us turn to the local static susceptibilities (||). In Fig. 6 the temperature dependences of Xn at SI = are 
displayed. Such dependences arc different at various sites. For example, the spin at site 2 of the chain shown in Fig. 
lb at low temperature shows little response to the applied field in contrast to spins at sites 1 and 3 (see Fig. 6b). 
It can be also noted that some spins in a nonuniform chain (e.g., the spins at sites 5-8 for the chain shown in Fig. 
lh) may exhibit almost a temperature-independent static susceptibility. Usually for the gapped (gapless) at SI = 
chains shown in Figs, lc, Id, le, If, lg (Figs, la, lb, lh) the low-temperature behavior of Xn and x differs only in 
quantitative details. However, in special cases we can observe a qualitative difference. Thus, the temperature behavior 
of the static susceptibility for a gapless chain mostly increases from a finite value to a maximum and then decreases 
inversely proportionally to temperature, whereas, for example, the local susceptibility at site 1 (2) of the chain shown 
in Fig. la (lb) exhibits a different behavior: as the temperature increases it decreases achieving a minimum at a finite 
temperature and then increases approaching the high-temperature asymptotic. 
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Figure 5: Local magnetizations m„ vs. field f2 at zero temperature for the chain shown in Fig. la for 5 = 0.3 (a), 0.6 (b), 0.9 
(c), 0.999 (d) (/ = 1). Solid curves correspond to the exact results, dotted curves correspond to the approximate ones obtained 
on the basis of Eqs. (jxs|) , (Q. 
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Figure 6: Local static susceptibilities \n at f2 = vs. temperature for the chains shown in Fig. 1 (7 = 1, 5 = 0.6). Panel a 



corresponds to a chain shown in Fig. la, panel b corresponds to a chain shown in Fig. lb, and so on. 
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6 Summary 



To conclude, we studied the local magnetizations and the local static susceptibilities of the spin-i XX chain in a 
transverse field with regularly alternating exchange interactions. These quantities can be calculated exactly using the 
Jordan- Wigner fermionization, Green function approach and continued fractions. We showed a relation between the 
alteration of the exchange interactions and i) the zero temperature local magnetizations along the chain at different ex- 
ternal fields and ii) the local susceptibilities along the chain at different temperatures. We found that the characteristic 
fields at which the zero temperature magnetization plateaus start and end up are the same for all sites contrary to the 
heights of plateaus, which are not universal but site-dependent and which depend on details of intersite interactions 
and the applied field. We interpreted the observed magnetization profiles from a viewpoint of the strong-coupling 
approach demonstrating a region of validity of that approximation. We discussed the temperature behavior of the 
local susceptibilities. 

The present study was supported by the DFG (projects 436 UKR 17/2/00 and Ri 615/6-1). O. D. acknowledges 
the kind hospitality of the Magdeburg University in the summer of 2000 when this paper was completed. 
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